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LONGITUDINAL  DYNAMIC  RESPONSES 

OF  A  STABILIZED  AIRPLANE 

By  Louis  H.  Smaus,  Marvin  R.  Gore, 
and  Merle  G.  Waugh 


SUMMARY 


The  dynamic  longitudinal  stability  of  an  airplane  with  autopilot 
was  predicted  by  combining  the  transfer  functions  of  the  autopilot 
obtained  from  tests  made  on  the  ground  with  those  of  the  airplane  meas¬ 
ured  in  flight  to  obtain  the  open-  and  closed-loop  frequency  responses 
and  transient  responses  for  the  combination.  These  predicted  responses 
were  then  compared  with  measured  flight  frequency  and  transient 
responses  for  three  airspeeds  and  various  autopilot  settings  of  dis¬ 
placement  and  rate  of  displacement  feedback.  The  analysis  procedures 
were  based  upon  linear  methods  and  agreement  was  good  when  elements 
were  operated  in  flight  within  the  linear  range  except  for  certain  con¬ 
ditions  in  which  rate  of  displacement  feedback  was  used. 


IRTRODUCTIOR 


The  purpose  of  the  investigation  was  to  determine  how  well  the 
longitudinal  dynamic  stability  of  an  autopilot-aircraft  combination 
could  be  predicted  from  the  separately  measured  characteristics  of  the 
autopilot  and  of  the  aircraft.  The  methods  of  analysis  were  based  on 
standard  servomechanism  theory  as  exemplified  by  reference  1.  A  gen¬ 
eral  survey  of  the  methods  of  predicting  the  dynamic  response  for  an 
autopilot-aircraft  combination  is  given  in  reference  2.  Most  of  the 
specific  equations  used  in  this  analysis  and  the  manner  of  diagramming 
the  closed-loop  system  were  developed  in  reference  3* 

The  transfer  functions  obtained  experimentally  for  the  aircraft 
and  for  the  autopilot  were  multiplied  together  and  the  resulting  open- 
loop  frequency  response  was  plotted  on  a  conventional  Nyquist  diagram 
to  indicate  the  relative  stability.  The  closed-loop  frequency  response, 
that  is,  the  ratio  of  pitch  response  to  a  sinusoidal  disturbance  in 
pitch  for  the  autopilot -aircraft  combination,  was  then  calculated 
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directly  from  the  open-loop  response.  In  addition,  the  transient 
responses  to  step  input  disturbances  were  calculated  by  an  approximation 
method.  These  three  predicted  responses  were  compared  with  those  meas¬ 
ured  in  flight ,  the  transient  response  and  the  closed-loop  frequency 
response  having  been  measured  directly  for  step  and  sinusoidal  inputs, 
respectively,  and  the  open-loop  response  derived  analytically  from  the 
closed-loop  response.  Calculations  and  measurements  were  made  for  air¬ 
speeds  of  85,  130,  and  200  knots  and  for  various  autopilot  settings  of 
displacement  and  rate  of  displacement  feedback. 


DEFINITIONS  AND  SYMBOLS 


Frequency  response:  A  frequency-dependent  vector  response  of  the  output 
of  a  system  to  a  sinusoidally  varying  input  function,  expressed  quanti¬ 
tatively  by  a  plot  of  amplitude  ratio  and  phase  angle  versus  frequency. 

Amplitude  ratio:  The  ratio  of  the  output  amplitude  to  the  input  ampli¬ 
tude.  For  a  closed-loop  system  this  is  ordinarily  nondimensiomalized 
by  dividing  by  the  amplitude  ratio  at  zero  frequency. 

Phase  angle:  The  angle  between  an  output  vector  and  input  vector.  When 
the  output  leads  the  input,  the  angle  is  positive. 

Transfer  function:  The  expression  defining  the  ratio  of  the  output  of 
a  component  to  the  input,  usually  expressed  as  a  complex  function  of 
the  frequency  variable  f . 

Closed-loop  response:  The  frequency  response  of  a  closed- loop  system, 
that  is,  one  which  possesses  feedback  and  is  sensitive  to  the  differ¬ 
ence  between  output  and  input. 

Open-loop  response:  The  frequency  response  of  an  open-loop  system. 

Servo  system:  That  part  of  the  autopilot  composed  of  the  amplifier  and 
servo  actuator  or  motor  and  its  internal  feedback  loop. 

Autopilot:  The  airplane  stabilizing  system  composed  of  the  servo  system 
and  the  feedback  gyros. 

Voltages,  angular  displacements,  and  transfer  functions  are  vector 
quantities  having  amplitudes  and  phase  angles  unless  otherwise  noted. 

aL  open-loop  transfer  function  of  autopilot -aircraft  combination 
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Ap  transfer  function  of  servo  system,  nondimens ionali zed 
Ar  transfer  function  of  rate  gyro 


transfer  function  of  aircraft  in  pitch 


2.718... 


frequency,  cycles  per  second 


follow-up  pickoff  constant,  volts  per  degree 

displacement  gyro  constant,  volts  per  degree 

static  control  gearing,  ratio  of  control  surface  deflection  to 
angular  displacement  input  to  autopilot,  degrees  per  degree 

rate  gyro  constant,  volts  per  cycle  per  second  per  degree 
oscillation  “ 


gain  of  follow-up  attenuator,  also  referred  to  as  sensitivity, 
percent 

gain  of  rate  gyro  attenuator,  percent 

amplitude  ratio  of  servo- system  frequency  response  AD, 
dimensionless  F 

amplitude  ratio  of  autopilot  frequency  response  when  rate  of  dis¬ 
placement  input  signal  is  included  (  Ap  — .PrAr  Y  dimensionless 

\  *g  / 

error  signal  of  servo  system  (input  to  amplifier)  when  the  servo 
system  is  tested  with  displacement  input  signal  only,  volts 

error  signal  of  servo  system  (input  to  amplifier)  when  the  servo 
system  is  tested  with  both  displacement  and  rate  of  displacement 
input  signals,  volts 

error  signal  of  servo  system  (input  to  amplifier)  when  the 
autopilot-aircraft  closed-loop  combination  is  tested,  volts 

error  signal  of  autopilot-aircraft  combination  (vj-vg),  volts 

feedback  voltage  of  servo  system,  volts 
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feedback  voltage  of  servo  system  when  tested  with  both  displacement 
and  rate  of  displacement  input  signals,  volts 

displacement  gyro  output,  volts 

input  signal  to  the  servo  system,  volts 

input  signal  to  the  autopilot -aircraft  combination,  volts 

rate  gyro  output,  modified  by  rate  attenuator,  volts 

control  surface  deflection,  degrees 

servo  displacement,  inches 

phase  angle  of  ve  relative  to  v^,  degrees 
phase  angle  of  ver  relative  to  v-^,  degrees 

phase  angle  of  servo-system  frequency  response  Ap,  vf  relative 
to  Vj_,  degrees 


phase  angle  of  Vf  relative  to  ve,  degrees - 
phase  angle  of  autopilot  frequency  response  ( Ap 


l+P  t»A 


when  rate  of  displacement  input  signal  is  included,  equivalent 
to  phase  angle  of  Vfr  relative  to  0,  degrees 

phase  angle  of  open-loop  autopilot-aircraft  combination  Al, 
(vg+vr)  relative  to  Vi  (and  (0+0r)  relative  to  0-jJ,  degrees 

phase  angle  of  vr  relative  to  0,  degrees 

angular  displacement,  attitude  of  aircraft,  degrees 

error  angle,  degrees  . 


hypothetical  input  angle  to  servo  system,  degrees 

input  angle  to  autopilot -aircraft  combination, 
degrees 

hypothetical  rate  feedback  angle,  degrees 


See  diagram 
page  11 
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DESCRIPTION  OF  TEST  EQUIPMENT 
Airplane 


The  airplane  used  in  these  tests  was  a  conventional  propeller- 
driven  Navy  dive  bomber  which  was  equipped  for  automatic  control.  A 
photograph  is  presented  in  figure  1. 

A  control  cable  and  pulley  system  connected  the  elevator  and  control 
stick.  The  autopilot  servo  actuator  was  attached  at  the  control  stick  end 
of  this  linkage. 


Autopilot 


Automatic  control  of  the  airplane  about  all  axes  was  furnished  by 
a  commercially  manufactured  autopilot  modified  to  include  rate  gyros. 

A  photograph  of  the  basic  components  of  the  pitch  channel  of  the  system 
is  shown  in  figure  2.  A  displacement  or  vertical  gyro  is  used  to  sense 
pitch  angle.  A  rate  gyro  senses  pitching  angular  velocity.  The  servo 
system  produces  an  elevator  deflection  in  accordance  with  the  gyro  out¬ 
puts  and  includes  the  components  described  as  follows.  An  amplifier 
converts  the  small  electrical  signals  received  from  the  gyros  and  follow¬ 
up  pickoff  to  currents  sufficiently  large  to  operate  the  solenoid- 
controlled  hydraulic  transfer  valve.  The  solenoid  transforms  the  elec¬ 
tric  current  to  a  mechanical  motion,  closing  or  opening  ports,  and  thus 
controlling  the  flow  of  hydraulic  fluid.  A  piston-type  servo  actuator 
converts  the  hydraulic  flow  to  a  linear  motion  having  sufficient  force 
to  actuate  the  control  surface  of  the  airplane.  The  actuator  is  con¬ 
nected  to  the  elevator-  cross  member  of  the  control-stick  mechanism.  A 
follow-up  pickoff  is  attached  to  the  piston  output  to  produce  an  elec¬ 
trical  signal  proportional  to  displacement  which  is  then  fed  back  to 
the  amplifier  to  complete  the  inner  servo  loop.  This  feedback  is  varied 
by  means  of  the  sensitivity  control,  a  potentiometer  controlling  the 
input  excitation  to  the  follow-up  pickoff.  Variation  of  this  control 
has  the  dual  effect  of  changing  the  dynamic  response  of  the  servo  system 
and  changing  the  ratio  of  surface  deflection  to  input  signal  which 
alters  the  response  of  the  autopilot -aircraft  combination. 

A  block  diagram  of  the  complete  autopilot -aircraft  loop  is  shown 
in  figure  3*  A  potentiometer  associated  with  the  rate  gyro  allows 
variation  of  the  amount  of  pitch-rate  signal  fed  back.  No  such  control 
is  provided  for  the  displacement  gyro  since  a  variation  of  the  servo 
follow-up  potentiometer  effectively  alters  the  pitch-angle  feedback. 
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Instrumentation 


^8ySoem  f°r  flisht  TOB  centered  about  a  six-channel 

by^Ls  of  roStSi. >erV°,.and  6leVator  surface  positions  were  measured 

^SL«es  ?itS  «nT  °rmerS  US6d  aS  electrical  Pick-offs  and  lever 

arm  xinicages .  Pitch  angle  was  measured  with  a  Sperry  A-12  vertical  j?vrn 

which  has  a  similar  pickoff .  The  alternating-current  voltages  from  thT 

j£o-m°n  “J?  at*i^de  pickoff s  were  demodulated  and  recorded  on  the 

oscillograph.  Rate  of  pitch  was  obtained  from  a  rate  gyro  with  a  micro 

dSectlv°o?  th^n400^ 7Cle  alternating-current  voltage  was  recorded 

the  oscillograph  without  rectification.  A  separate  galva¬ 
nometer  was  used  to  record  the  current  to  the  solenoid  valve,  this  beimr 

^  TOltage  to  amplifier.  A  ataXl  Mof  ^ 
airspeed-altitude  recorder  also  was  used. 


For  sinusoidal  response  tests  the  sine-wave  input  signal  was 
obtained  from  a  device  in  which  a  rotating  selsyn  was  driven  by  a  ball- 
disc  variable-speed  mechanism.  The  amplitude  of  the  kOO-cycle  output 
voltage  from  the  selsyn  could  be  varied  in  frequency  from  0  to  10  cycles 
per  second.  A  contactor  was  operated  once  per  cycle  at  a  zero  output 
voltage  point,  introducing  a  signal  to  the  oscillograph  which  provided 

a  zero-phase  reference  marker.  The  sine-wave  amplitude  was  set  and 
measured  on  the  ground. 

For  transient-response  tests  pulses  were  introduced  by  a  motor- 
driven  potentiometer  arrangement.  Both  amplitude  and  time  base  could 
be  adjusted  on  the  ground.  Steps  were  obtained  by  a  simple  switch.  The 
various  input  voltages  were  introduced  in  series  with  the  other  elements 

in  the  autopilot  signal  circuit  and  are  indicated  by  the  symbol  vT  in 
rigure  3«  - 

,  ^°r  ground  ^sts  of  the  servo  system  a  Brush  recorder  was  used 

to  indicate  servo  and  surface  positions,  and  the  zero-phase  reference 
marker  was  superimposed  on  one  of  the  traces. 


THEORETICAL  ANALYSIS 


method  of  predicting  the  dynamic  response  of  the  stabilized 
(autopiiot-controHed)  airplane  is  based  on  the  determination  and  combi¬ 
nation  of  the  transfer  functions  of  the  aircraft  and  of  the  autopilot. 
The  relative  stability  may  then  be  ascertained  from  an  inspection  of  the 
resultant  open-ioop  frequency  response,  closed-loop  frequency  response, 
and  transient  response.  These  responses  are  mathematically  related  as 

fabse<luently>  ^t  each  provides  information  on  the  relative 
stability  not  readily  apparent  from  the  others. 
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The  Aircraft 


Transient  experimental  flight  data  was  the  most  important  source 
of  aircraft  frequency-response  characteristics.  This  information  was 
obtained  from  flight  film  records  by  a  method  of  analysis  based  upon  the 
Fourier  transform  as  outlined  as  follows. 

For  many  considerations  of  dynamic  stability  the  longitudinal 
motions  of  an  airplane  can  be  closely  described  by  a  second-order  differ¬ 
ential  equation  (reference  ^) .  The  quantities  involved  in  the  stabili¬ 
zation  system  of  the  test  airplane,  namely,  the  angle  of  pitch,  Q,  as 
controlled  by  the  elevator  position,  6,  may  be  related  by  the  equation 


d 29 
dt2 


+  b 


—  +  k9  =  Cxb  +  C0 
dt 


5  dt 


This  differential  equation  may  be  converted  by  the  Laplace  transfor¬ 
mation  into  its  transfer-function  form: 


e(s)  ,_Ji£i22_s(s) 

s(  s2+bs+k) 

When  6  and  5  are  known  experimentally  as  functions  of  time, 
they  can  be  converted  into  the  s  plane  by  the  Laplace  integrals 


(1) 


0(s)  =  r“  0(t)e“stat 
5(s)  =  Pm  8(t)e"stat 

where  0(t)  and  6(t)  are  assumed  to  be  zero  for  all  t<0. 
Equation  (l)  can  then  be  written 

fW  9(t)e-stat  =  ._Cib+Co  r“8(t)e-stat 

^  s(s2+bs+k)  J 
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The  frequency  response  may  now  he  calculated  by  replacing  s  by  jw, 
where  to  is  any  arbitrary  value  of  angular  frequency.  (See  reference  1.) 
Thus,  frequency  response  is  written: 


0(  dw) 

_  jo))+C0 

-P 

3 

**“D 

1 

CD 

, 

CD 

8 

6(  Jw) 

( «)“)[■'(  jw)S+b(  jto)+k] 

A 

The  two  integrals  are  called  Fourier  integrals  and  are  usually 
directly  calculable.  When  these  integrals  do  not  converge,  limiting 
values  may  be  obtained  by  replacing  the  integrals  with  the  following  more 
general  form  of  the  Laplace  transformation  in  which  s  =  0  +  jto: 

0(ju)  =  CT^mo  y'<X,0(t)e-(o+j“)tdt 

O 

Thus  by  equation  (2),  when  Q  and  5  are  known  as  functions  of 
time,  the  frequency  response  may  be  computed  for  any  arbitrary  value  of 

angular  frequency,  w,  by  evaluating  the  two  simple  integrals  above. 

*  ' 

Considerable  work  can  be  saved  if  0(t)  and  5(t)  are  chosen  so 
that  the  transient  portions  are  short.  (See  fig.  4.)  The  integrals 
then  involve  statistical  integration  up  to  some  time,  T,  when  steady 
state  is  reached  and  an  analytical  expression  will  finish  the  evaluation 
of  the  integral  from  T  to  infinity.  The  statistical  work  is  begun  by 
dividing  the  interval  t  =  0  to  t  =  T  into  increments  of  about  0.1 
second,  and  the  integration  is  performed  numerically,  utilizing  some 
approximation  system  such  as  Simpson's  Rule. 

This  method  for  obtaining  frequency-response  functions  from 
transient  data  is  discussed  in  greater  detail  in  reference  5* 


The  Autopilot 


The  characteristics  of  the  autopilot  and  its  components  were 
obtained  primarily  by  tests  on  the  ground  of  the  equipment  while 
installed  in  the  airplane.  Since  the  predicted  autopilot -aircraft 
response  is  theoretically  valid  only  if  the  individual  component 
responses  are  linear,  the  autopilot  was  examined  with  this  consideration 
in  mind. 
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The  servo  system  is  linear  in 
operation  only  for  input  signals  to 
the  amplifier  not  exceeding  a  cer¬ 
tain  value.  This  nonlinearity 
results  from  saturation  of  the 
amplifier  which  has  a  static  input- 
output  relation  as  shown  in  the 
sketch.  The  amplifier  input  voltage 
is  more  coveniently  referred  to  as 
the  servo -system  error  voltage,  ve, 
as  can  be  seen  from  figure  3«  Thus 
it  becomes  necessary  to  know  the 
error  voltage  to  determine  whether 
the  servo  system  is  operating  within 
its  linear  range.  The  error  voltage 


input  voltage 


is,  of  course,  a  function  of  the  input  voltage  to  the  servo  system  v± 
as  well  as  of  the  feedback  voltage  vp  from  the  servo  output. 


If  the  response  of  the  servo  system  alone  (without  the  gyros)  is 
being  considered,  the  input  to  the  system  is  v-j_.  An  expression  for  the 
error  voltage  in  this  case  is  given  in  reference  3,  equation  (5),  and  is 


where  R  and  £f  are,  respectively,  the  nondimens ional  amplitude  ratio 
and  phase  angle  of  the  closed-loop  servo-system  response  Ap. 

When  the  response  of  the  autopilot,  that  is,  the  control- surface 
response  6e  to  an  angular  attitude  input  G  to  the  gyros,  is  con¬ 
sidered,  the  error  voltage  ver  in  this  case  is,  from  equation  (12)  of 
reference  3> 


where  Rfr  and  €fr  are,  respectively,  the  nondimens ional  amplitude 
ratio  and  phase  angle  of  the  autopilot  response  when  both  rate  and  dis¬ 
placement  gyros  are  oscillated.  The  term  Prkpf  is  the  amplitude  of 
the  transfer  function  Pr-A-r  for  the  rate  gyro.  The  displacement  gyro 
output,  for  practical  purposes,  is  related  to  the  input-  Q  by  a 
constant  kg. 


From  a  knowledge  of  the  servo-error  voltage  it  is  possible,  in  many 
instances,  to  choose  input  magnitudes  small  enough  to  insure  operation 
of  the  servo  in _ the  linear  range.  This  is  not  always  the  case,  however. 


10 


KA.CA  TF  2578 


and  several  of  the  test  results  presented  herein  involve  the  nonlinear 
range,  but  no  attempt  is  made  to  calculate  the  effects  of  nonlinearity 
on  the  predicted  results. 

A  second  source  of  nonlinearity  exists  in  the  elevator  linkage. 
Coulomb  friction  in  the  control-surface  hinge  and  flexing  of  the  pulley 
brackets  and  supporting  deck  are  the  probable  causes. 


Stability  Criteria 


Several  criteria  exist  for  determining  stability.  One  of  the  most 
popular  is  based  on  the  use  of  the  Nyquist  diagram  which  is  explained 
fully  in  reference  1  and  other  references  listed  therein.  Briefly,  a 
polar  plot  in  the  complex  plane  is  made  of  the  open-loop  frequency  x 

response.  For  the  cases  considered  in  this  investigation,  an  encircle-  >5 

ment  of  the  -1  4-  JO  point  represents  an  unstable  system,  whereas  coin¬ 
cidence  with  the  point  represents  a  condition  of  neutral  stability.  The 
latter  case  is  simply  equivalent  to  saying  that  the  output  magnitude  is 
the  same  as  the  input  (the  system  gain  is  unity),  but  is  l80°  out  of 
phase.  Since  the  displacement  feedback  in  a  servo  system  is  negative  * 

and,  hence,  has  an  additional  l80°  phase  lag,  it  will  lag  a  total  of  36O0 
and  therefore  add  to  the  input  signal.  Therefore,  any  oscillations 
which  may  start  are  self-sustaining  and  the  system  will  hunt  indefinitely  *>' 

and  with  constant  amplitude . 

The  nearness  of  the  plot  to  the  -1  4  jO  point  is  an  indication  of 
the  relative  stability  of  the  system  and  is  often  given  in  terms  of  gain 
and  phase  margins.  Phase  margin  is  defined  as  the  angle  between  l80° 
and  the  point  at  which  the  open-loop  response  passes  through  the  circle 
of  unit  magnitude.  Gain  margin  is  the  reciprocal  of  the  open-loop 
response  magnitude  when  the  phase  angle  is  l80°.  (See  reference  1.) 

However,  no  simple  quantitative  relations  correlate  the  Nyquist  plot 
with  characteristics  of  the  transient  response  such  as  period  and  time 
to  damp  to  a  specified  amplitude,  except  for  a  second-order  system. 

The  closed-loop  frequency  response  plot  also  gives  an  indication  of 
the  stability.  The  amplitude  of  the  resonant  peak  is  a  rough  measure  of 
the  damping,  and  the  resonant  frequency,  if  it  clearly  exists,  is  close 
to  the  transient  frequency  and  the  undamped  natural  frequency  of  the 
system. 

Perhaps  the  final  measure  of  stability  is  the  transient  response 
itself.  Usually  a  particular  form  of  this  response  is  desired,  often 
one  that  has  an  equivalent  damping  ratio  in  the  neighborhood  of  0.6 
cr^"^^cal  damping  and  a  certain  speed  of  response  or  natural  frequency. 

These  criteria  may  be  modified  by  other  limitations,  such  as  the 
maximum  allowable  acceleration  that  may  be  imposed  on  personnel  or 
airframe . 


/ 


/ 


NACA  TN  2578 


11 


The  Autopilot -Aircraft  Loop 


A  block  diagram  of  the  autopilot -aircraft  loop  is  given  in  figure  3* 
For  purposes  of  analyzing  the  open-loop  characteristics,  the  loop  may  be 
opened  anywhere,  but  it  is  most  convenient  to  break  it  between  the  air¬ 
craft  output  and  the  gyro  input.  The  open-loop  transfer  function  is 
then  given  in  terms  of  a  ratio  of  the  pitch-angle  output  to  a  pitch-angle 
input.  This  is  the  function  that  is  plotted  on  a  Nyquist  diagram.  As 
previously  stated,  varying  the  amount  of  elevator  deflection  for  a  given 
vertical  gyro  displacement  changes  the  gain  of  the  over-all  system.  The 
change  in  amplitude  on  the  Nyquist  plot  is,  of  course,  proportional  to 
the  change  in  system  gain,  and,  hence,  directly  affects  the  stability  of 
the  closed-loop  combination. 

A  closed-loop  response  is  obtained  by  inserting  an  input  signal  at 
some  point  and  measuring  an  output  response  at  any  other  point  with  the 
loop  closed.  The  output  generally  of  interest  is,  again,  the  pitch 
angle.  It  is  then  desirable  to  use. a  pitch  input,  but  in  flight  it  is 
impractical  to  feed  in  a  sinusoidally  varying  pitch  angle.  Therefore, 
an  equivalent  pitch  input  is  obtained  by  inserting  a  voltage,  Vj,  in  the 
autopilot  signal  circuit  (fig.  3)*  This  voltage  is  related  to  the  hypo¬ 
thetical  pitch  input  then  by  the  same  constant,  kg,  relating  voltage 
output  per  degree  input  for  the  displacement  gyro.  The  input  for  a 
transient  response  is  introduced  in  the  same  manner. 


Fundamental  Eelations 


Relations  between  open-loop  and  closed-loop  frequency  responses  are 
given  generally  in  servomechanisms  texts.  The  particular  forms  used 
herein  were  derived  in  reference  3  for  application  to  the  autopilot- 
aircraft  combination.  They  are  expressed  in  terms  of  the  quantities 
actually  measured  during  tests.  In  this  connection,  as  is  shown  in 
reference  3>  it  is  helpful  to  redraw  figure  3  as  follows: 
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The  equivalent  angular  input  is  (where  vj  4  kg  0j)  and  the  error 
angle  is  The  term  kp  represents  the  autopilot  gearing,  the  static 

ratio  of  elevator  deflection  to  vertical  gyro  displacement.  The  nondi- 
mensional  frequency  response  of  the  autopilot  with  the  output  measured 
at  the  control  surface  and  without  pitch-rate  feedback  is  Ap.  The 
relative  amount  of  rate  to  displacement  feedback  is  represented  by  the 
term  PrAr/kg  where  Ar  is  the  rate -gyro  transfer  function,  Pr  the 
rate  attenuation  factor,  and  kg  the  vertical  gyro  constant.  The  air¬ 
craft  response  is  A q. 

It  should  be  noted  that  all  quantities  are  considered  as  vectors, 
possessing  both  magnitude  and  phase  angle,  unless  otherwise  noted.  The 
term  vj,  or  its  equivalent  0j,  is  the  reference  and  has  zero  phase 
angle . 

The  equations  presented  in  the  following  paragraphs  were  used  in 
the  calculations  involved  in  predicting  and  analyzing  the  performance 
of  the  autopilot -aircraft  combination.  They  are  presented  without 
formal  derivation  and  may  be  derived  from  the  preceding  diagram  and 
elementary  servo  theory.  For  further  details  of  these  relations  and 
others  governing  the  calculation  of  servo-error  voltage,  servo  response 
for  a  change  of  servo  gain,  and  servo  response  with  addition  of  rate  of 
displacement  input  signal  the  reader  is  referred  to  reference  3- 

Open-loop  frequency  response  from  closed-loop  frequency  response. - 
The  open-loop  response  Ap  may  be  calculated  from  the  closed-loop 
response  0/6 j  as  measured  in  flight  from  the  equation 


From  this  it  can  be  seen  that  in  addition  to  measuring  the  closed- 
loop  response  it  is  necessary  to  evaluate  the  feed-back  factor 
(l  +  PrAr/krr)  .  This  can  be  done  in  two  ways.  If  the  transfer  functions 
kg  and  PrAr  for  both  displacement  and  rate  gyros,  respectively,  are 
known,  the  factor  may  readily  be  calculated  for  the  frequency  range  of 
interest. 

On  the  other  hand,  if  the  frequency  response  of  the  autopilot  servo 
system  is  known  for  a  rate  setting  of  zero  and  for  the  rate  signal  being 
considered,  the  feed-back  factor  may  be  obtained  from  the  two  responses. 
It  may  readily  be  seen  from  the  foregoing  sketch  that  the  nondimens ional 
autopilot  dynamic  response  with  rate  signal  present  is  the  product  of 
the  servo-system  response  Ap  and  the  feed-back  factor  (l  +  PrAr/kg). 
(The  gearing  term  kp  shown  in  the  sketch  is  simply  a  constant  which 
converts  the  nondimensional  response  to  its  absolute  magnitude.) 
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Therefore  the  feed-back  factor  can  be  evaluated  from  the  ratio  of  the 
rate  response  to  no-rate  response,  or 


measured 


(6) 


Open-loop  frequency  response  from  autopilot  and  aircraft  component 
responses,-  The  open-loop  response  is  simply  the  "product  of  the  indi- 
vidual  transfer  functions  around  the  loop. 


AL  =  kp  Ap  ^1  +  A e  ( 7) 

As  was  shown  previously,  the  term  Ap(l  +  PrAr/kg)  represents  the 
nondimens ional  response  of  the  autopilot  servo  system  with  raue-of- 
displacement  signal  added  and  may  be  measured  directly.  Or,  if  more 
convenient,  the  feed-back  factor  (l  +  PrAr/kg)  may  be  calculated  from 
the  gyro  transfer  functions. 

Closed-loop  frequency  response  from  open-loop  frequency  response. - 
The  closed-loop  response  0/0  j  cannot  be  calculated  directly  from  the 
open-loop  response  alone  but  must  take  into  account  the  feedback  due  to 
the  rate  gyro.  Hence, 

0  _  kpApAp _ 

01  1+kpApAg  ^1+ 

§ 

In  terms  of  the  open-loop  response  from  equation  (7)  and  the  feed-back 
factor,  the  above  equation  may  be  rewritten  as 


(9) 


Error  voltage  for  the  autopilot -aircraft  combination.-  The  servo - 
system  error  voltage  has  been  discussed  previously  in  connection  with 
tests  of  the  components .  When  it  is  desired  to  predict  the  autopilot- 
aircraft-combination  response,  it  is  again  necessary  to  calculate 
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the  error  voltage  for  the  input  signal  being  considered  in  order  to 
determine  whether  the  amplifier  is  operating  within  its  linear  range. 

The  expression  for  the  error  voltage  in  this  case,  designated  vec,  is, 
from  reference  3,  equation  (42), 


vec I  =  VI 


Transient  response  from  the  closed-loop  frequency  response.-  The 
transient  response  for  a  step-input  disturbance  was  obtained  from  the 
closed-loop  response  by  Floyd's  Method,  an  approximation  method  which 
is  explained  in  detail  in  chapter  11,  reference  1. 

Briefly,  the  inverse  transform  h(t)  of  H(s),  where  h(t)  is  the 
transient  response  to  an  impulse  and  H(s)  is  a  function  of  the  complex 
operator  s,  is 


h(t)  =  S  ^s^etSds  (ll) 

C-joo 

Under  the  conditions  where  s  may  be  replaced  by  ju,  an  exact  expres¬ 
sion  for  the  inpulse  response  in  terms  of  the  real  part  of  the  closed- 
loop  frequency  response  H( ju)  is 


h(t)  =  2  J  [Re  H(jw)  cos  tu]  du  (12) 

o 

This  integral  may  be  approximated  graphically.  The  procedure  is  to 
plot  Re  H(jto)  against  cj  and  approximate  the  exact  shape  of  the  curve 
by  a  series  of  straight-line  segments.  This  straight-line  approximation 
defines  a  series  of  trapezoidal  functions  each  of  which  can  be  evaluated 
by  equation  (12).  The  approximate  value  of  h(t)  is  then  obtained  by 
adding  the  resultant  time  functions. 
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For  a  series  of  trapezoidal  functions  such  as  is  shown. 


Re  H(jw) 


equation  (12)  can  be  evaluated  as 


K 


n 


k 

An 

u)n 


number  of  trapezoids  used  in  approximating  curve 
area  of  n^*1  trapezoid  (rncon) 

2 


An 


2 


The  step  response  may  then  he  obtained  by  a  graphical  integration 
with  respect  to  time  of  the  impulse  response. 


TEST  RESULTS  AND  CALCULATIONS 


In  this  section  are  presented  the  test  data  and  associated  calcula¬ 
tions  required  for  analysis  of  system  stability  at  the  three  test  air¬ 
speeds-  of  85,  130,  and  200  knots  indicated  at  an  altitude  of  10,000  feet. 
The  aircraft,  servo  system,  and  autopilot  gearing  are  first  considered 
individually.  Hinge  moment,  required  for  gearing  calculations,  is  also 
considered.  The  flight  closed-loop  frequency  responses  and  their  corre¬ 
sponding  open-loop  responses  are  then  presented  and  compared  with  open- 
and  closed-loop  responses  calculated  from  the  component  transfer 
functions.  Finally,  transient  responses  predicted  from  component 
frequency  responses  are  compared  with  those  measured  in  flight  for  the 
closed-loop  system. 
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Aircraft  Transfer  Function 


The  dynamic  characteristics  of  this  component  of  the  loop  were 
obtained  by  analysis  of  flight-test  data.  Two  types  of  flight  maneuver, 

wave  forced  oscillation  and  transient,  were  used  as  the  source  of 
aircraft-frequency-response  data.  The  inputs  were  voltage  disturbances 
of  prescribed  form  introduced  into  the  autopilot -aircraft  loop  (indi¬ 
cated  by  vj  in  fig.  3).  This  technique  avoided  the  considerable  effort 
necessary  to  perform  such  forced  oscillations  in  an  unstabilized  aircraft 
necessitating,  here,  addition  of  a  simple  sine-wave  voltage  generator 
only.  Transient  inputs  introduced  in  this  fashion  may  be  considered 
unusual  to  those  who  have  not  worked  with  stabilized  aircraft.  The 
actual  elevator  movement  that  produces  the  aircraft  response  does  not 
assume  the  step  function  or  pulse-type  form,  but  is  itself  a  damped 
oscillation.  Thus,  the  work  involved  in  analyzing  the  records  is 
increased.  The  quantities  measured  for  the  study  of  this  aircraft,  which 
is  stabilized  in  pitch,  were  pitch  angle  6  and  elevator  input 
angle  6e .  Sample  forced-oscillation  and  transient-response  records  are 
presented  in  figure  4. 

A  straightforward  method  was  utilized  in  extracting  the  aircraft 
frequency-response  characteristics  from  the  flight  film  records  for  the 
sine-wave  oscillation  tests.  While  the  9  and  6e  records  are  not 
true  sine  waves,  they  were  considered  near  enough,  in  most  of  the  data, 
to  warrant  a  simple  analysis.  The  amplitude  ratio  was  obtained  by  meas¬ 
urement  of  the  amplitudes  of  the  peaks  and  phase  angle  from  the  time 
differences  between  corresponding  intercepts  of  the  line  drawn  to  equal¬ 
ize  the  half-cycle  time  intervals .  Records  not  lending  themselves  to 
this  procedure  represented  such  a  small  portion  of  the  data  that  they 
were  disregarded. 

The  frequency  response  was  obtained  from  transient  data  by  use  of 
equation  (2)  as  discussed  previously. 

The  final  frequency-response  characteristics  are  a  combination  of 
the  results  obtained  from  both  types  of  flight  tests.  Sine-wave- 
o^c illation  data  provide  amplitude  ratio  values  that  repeat  within 
±3  percent  for  low  values  of  frequency  (0-0.3  cps)  and  to  greater  accu- 
racy  for  values  up  to  1.6  cps.  Sine-wave-oscillation  phase  angles  were 
not  considered  to  be  of  usable  accuracy  for  the  types  of  analyses  that 
were  attempted.  Transient  analyses  provided  consistent  phase  angles 
and  amplitude  ratios  that  varied  within  a  maximum  deviation  of  about 
5  percent.  The  mean  values  of  the  transient-response  amplitude  ratios 
agreed  well  with  those  of  the  frequency  response  and  the  final  amplitude- 
ratio  curves  (fig.  5(&))  are  the  faired  average  of  both  sine-wave 
oscillations  and  transient  results.  The  final  phase-angle  character¬ 
istics  (fig.  5(b))  are  exclusively  from  transient  data. 
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As  a  result  of  this  investigation,  it  appears  that,  from  considera¬ 
tions  of  flight  time  and  practicable  analysis  procedures,  the  transient 
flight  test  is  superior  for  obtaining  frequency-response  data  for  the 
airplane  transfer  function. 


Autopilot  Frequency  Response 


Extensive  ground  tests  were  made  of  the  autopilot.  These  were 
performed  with  most  of  the  equipment  remaining  in  the  fuselage  to  simu¬ 
late  the  actual  flight  setup  as  closely  as  possible.  The  flight  aero¬ 
dynamic  loads  upon  the  elevator  were  simulated  by  torsional  springs 
producing  the  required  hinge  moments.  The  autopilot  frequency-response 
characteristics  were  determined  for  several  values  of  displacement  and 
rate -of -displacement  feedback  and  input-signal  amplitudes. 

In  order  to  choose  input  signals  of  a  magnitude  low  enough  to  insure 
linear  operation  of  the  servo  amplifier,  it  was  first  necessary  to  meas¬ 
ure  the  static  characteristic  of  the  amplifier.  This  relation  is  shown 
in  figure  6  where  the  output  current  is  plotted  as  a  function  of  the 
input  voltage  to  the  amplifier  (servo-system  error  voltage).  It  can  be 
seen  that  the  relation  is  linear  within  about  10  percent  over  a  range 
of  ±0*35  volt  about  the  voltage  value  required  for  zero  unbalance 
current.  This  value  of  0.35  volt  is  designated  the  nonlinearity  level 
although  it  is  evident  that  the  system  does  not  depart  rapidly  from 
linearity  for  another  0.1  volt  or  so. 

Dynamic  tests  were  conducted  on  the  autopilot  servo  system  coupled 
to  the  elevator  control  £urface  first  using  an  electrical  sinusoidal 
input  signal.  An  amplitude  ratio  expressed  as  the  ratio  of  the  output 
motion  at  any  frequency  to  that  at  zero  frequency  and  a  phase  angle 
representing  the  number  of  degrees  the  output  motion  leads  (considered 
plus)  or  lags  (minus)  the  input  signal  were  obtained.  Tests  were  made 
for  several  values  of  simulated  hinge  moment  but  it  was  found  that  the 
dynamic  response  did  not  differ  materially  over  the  frequency  range  of 
primary  interest,  0-1  cps.  Hence  the  no-load  responses  were  used  for 
the  analysis.-  Loading  of  the  control  surface  does  change  the  gearing, 
however,  and  this  is  discussed  in  the  following  two  sections.  The  no- 
load,  nondimens ional  frequency  reponse  Ap  for  the  servo  system  is 
presented  in  figure  7(a)  for  a  range  of  sensitivity  settings,  namely, 

24,  33*  42,  52,  and  63  percent.  The  magnitude  of  the  input  signal  was 
±0.115  volt,  corresponding  to  about  ±l/4°  in  pitch  for  the  gyro  constant 
used  in  the  analysis.  This  magnitude  was  low  enough  to  allow  linear 
operation  of  the  servo  system  throughout  the  frequency  range. 

To  determine  the  response  for  rate  signal  in  addition  to  displace¬ 
ment  signal,  tests  were  conducted  with  the  gyros  mounted  on  a  sinus¬ 
oidally  oscillating  table  and  their  combined  electrical  output  fed  to 
the  servo  system.  This  nondimens ional  response  is  designated  by  the 


18 


MCA  TN  2578 


factor  Ap(l+PrAr/kg)  and  is  shown  in  figure  7(h)  for  a  sensitivity  of 
2k  percent  and  rate  settings  of  0,  8,  20,  and  31  percent.  (These  rate 
settings  gave  values  for  the  ratio  of  rate  to  displacement  signals, 
PrAr/kg,  of  0.83f,  2.07f,  and  3 • 21f ,  respectively,  up  to  a  frequency 
of  1.2  cps.  At  higher  frequencies  the  amplitude  of  Ar  departed  from 
its  linear  relationship  with  frequency. )  The  input  magnitude  was  ±l/k° 
of  table  oscillation.  With  this  input  the  responses  with  rate  signal 
remained  linear  up  to  a  little  more  than  1  cps.  The  responses  for  all 
rate  signal  values  reached  the  nonlinearity  level  between  1  and  2  cps, 
the  response  at  the  lowest  rate  value  being  linear  almost  to  2  cps. 

The  response  with  rate  signal  was  also  calculated  from  the  measured 
values  of  Ap,  Pr,  Ar,  and  kg  for  the  same  conditions  as  above.  The 
agreement  with  the  response  obtained  from  the  oscillating  table  tests 
was  very  good.  (See  reference  3 • )  Hence,  for  a  sensitivity  of 
k2  percent,  the  rate  responses  were  calculated  rather  than  measured 
directly. 

.  Hinge  Moment 


Ground  tests  of  the  autopilot  installed  in  the  airplane  disclosed 
the  fact  that  hinge  moment  directly  affects  the  gearing  factor  kp. 

With  a  flexible  linkage  connecting  the  servo  actuator  and  control  sur¬ 
face,  as  indicated  in  figure  3 ,  the  gearing  will  be  decreased  by  addi¬ 
tion  of  any  load  on  the  elevator  due  to  stretching  of  the  control  cable. 
The  hinge  moment  was  determined  from  flight  tests  in  order  to  eliminate 
it  as  a  possible  source  of  error  in  the  predicted  autopilot -aircraft 
responses . 

The  spring  properties  of  the  linkage  between  servo  actuator  and 
control  surface  were  used  to  determine  the  hinge  moments  encountered 
in  flight  at  the  three  test  airspeeds.  The  linkage  was  calibrated  cn 
the  ground  and  its  spring  constant  determined.  Both  servo  and  surface 
positions  were  recorded  on  the  ground  under  no  load  and  in  flight  at 
the  three  airspeeds  during  the  course  of  the  frequency-response  tests. 
Thus  the  change  in  the  ratio  of  surface  to  servo  deflections  between 
ground  and  flight  conditions  was  a  measure  of  the  hinge  moment  under 
dynamic  tests.  The  average  value  of  a  number  of  runs  for  each  indicated 
airspeed  was  as  follows: 


Foot-pounds 

Knots 

per  degree 

85 

3.0 

*9 

130 

10.0 

200 

21.6 
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Hinge  moments  were  also  estimated  from  wind-tunnel  measurements  on 
a  similar  airplane  and  were  sufficiently  close  to  the  actual  values  as 
to  have  made  no  difference  in  the  determination  of  the  gearing  factors. 


Autopilot  Gearing 


The  autopilot  gearing  factor  kp  is  defined  as  the  static  ratio 
of  the  control-surface  deflection  to  pitch-attitude  input  to  the  auto¬ 
pilot.  In  this  particular  installation  it  is  controlled  by  the  sensi¬ 
tivity  potentiometer  but,  as  shown  previously,  it  is  also  a  function  of 
hinge  moment  and,  hence,  airspeed. 

The  gearing  factor  was  determined  in  several  ways,  all  of  which 
gave  somewhat  different  results.  The  most  complete  determination  was  by 
means  of  introducing  an  electrical  signal  of  several  magnitudes,  corre¬ 
sponding  to  various  pitch  angles,  to  the  servo  system  with  the  airplane 
on  the  ground  and  measuring  the  control-surface  deflection.  This  was 
done  for  no  load  on  the  surface  and  at  three  values  of  simulated  hinge - 
moment  load,  2,  8,  and  20  foot-pounds  per  degree,  corresponding  roughly 
to  the  test  airspeeds  of  85,  130,  and  200  knots.  A  range  of  sensitivity 
settings  was  also  covered  for  each  load.  To  obtain  the  corresponding 
pitch  angle,  the  vertical  gyro  was  calibrated  by  rotating  it  0.1°  at  a 
time  and  measuring  the  output  voltage.  Thus  the  gyro  constant  kg  was 
found  to  be  0.51  volt  per  degree.  The  gearing  is  then  the  product  of 
this  constant  and  the  values  of  surface  deflection  per  volt  input  to  the 
servo  system.  To  obtain  a  linear  relation,  the  reciprocal  of  the 
gearing  l/kp  is  plotted  in  figure  8  against  sensitivity  for  the  various 
loads.  The  data  are  replotted  in  figure  9  to  show  the  reciprocal 
gearing  l/kp  plotted  against  load  for  varying  sensitivities.  From 
this  figure  the  gearings  corresponding  to  the  test  airspeeds  were 
obtained  and  used  in  the  analysis  presented  in  this  report. 

The  gearing  for  the  no-load  case  was  determined  also  by  measuring 
the  control -surface  deflection  for  various  attitude  settings  of  the 
vertical  gyro  which  was  connected  to  the  rest  of  the  autopilot  in  the 
normal  manner  as  in  figure  3*  In  another  test  the  gyro  was  mounted  on 
an  oscillating  table  and  oscillated  sinusoidally  at  low  frequencies 
(between  0.1  and  0.2  cps)  at  several  amplitudes,  and  the  control- surface 
deflections  were  recorded.  The  values  obtained  from  these  two  types  of 
tests  were  not  in  good  agreement  with  each  other  or  those  obtained 
previously.  Furthermore,  it  was  observed  that  the  values  seemed  to 
depend  on  the  servo  displacement,  being  considerably  less  for  small  dis¬ 
placements.  This  is  the  type  of  behavior  typifying  systems  containing 
nonlinearity  caused  by  play  or  backlash,  and  indicates  one  probable 
cause  for  discrepancies  between  flight  and  predicted  responses. 
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This  effect  could  not  be  resolved  without  considerably  more 
investigation  than  was  justified  in  this  case.  The  values  obtained 
from  figure  9  gave  predicted  autopilot -aircraft  open-loop  responses 
.which,  on  the  whole,  most  nearly  matched  those  obtained  from  flight  and 
therefore  were  used  in  the  analysis. 


The  Autopilot-Aircraft  Combination 

Open-loop,  closed-loop,  and  transient  responses  for  the  autopilot- 
aircraft  combination  were  predicted  from  tests  of  the  components  and 
also  obtained  experimentally  in  flight  over  a  wide  range  of  sensitivity 
and  rate-of -pitch  feedback  conditions  for  each  of  the  three  airspeeds. 
The  conditions  tabulated  were  selected  for  presentation  in  this  report 
not  only  because  they  include  representative  sensitivity  and  rate  set¬ 
tings,  but  also  because  they  illustrate  the  effect  of  changes  in  air¬ 
speed,  sensitivity,  and  rate  upon  the  system  performance . 


Airspeed 

Sensitivity 

Rate 

Figure 

(knots) 

(percent) 

(percent) 

(number) 

85 

42 

0 

10(a) 

130 

42 

0 

10(b) 

200 

42 

0 

10(c) 

130 

63 

0 

10(d) 

130 

24 

8 

10(e) 

130 

24  , 

31 

10(f) 

130 

42 

8 

10(g) 

130 

42 

20 

10(h) 

130 

42 

31 

10  ( i) 

The  open-loop,  closed-loop,  and  transient  responses  each  provide 
information  of  a  vital,  although  different,  nature  about  the  stability 
and  performance  of  the  autopilot-aircraft  combination  and  is  plotted 
upon  the  same  page  for  each  of  the  above  conditions.  (See  figs.  10(a) 
to  10(i).)  In  each  figure  the  predicted  response  is  compared  with  the 
response  measured  in  flight  or  one  derived  therefrom.  The  methods  of 
calculation  have  been  presented  in  the  foregoing  sections  and  a  sample 
calculation  for  each  type  is  given  in  the  appendix. 

All  calculations  were  based  on  the  assumption  of  linear  operation 
aH  components.  This  assumption  holds  throughout  the  frequency  range 
for  very  small  inputs  only.  For  the  flight  tests  an  input  magnitude 
equivalent  to  ±1/2  was  used.  This  value  was  not  low  enough  to  insure 
linear  operation  all  the  time,  but  it  was  as  small  as  practicable  from 
the  standpoint  of  accuracy  of  measurements.  It  is  also  believed  to  be 
a  reasonable  value,  one  that  might  likely  be  encountered  in  flight  as 
an  external  disturbance.  The  error  voltage  for  the  closed-loop  combina¬ 
tion  was  calculated  in  each  case . 
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The  open-loop  frequency  response  is  shown  at  the  top  of  each  page 
plotted  on  polar  coordinates  representing  amplitude  and  phase  angle. 

This  form  is  commonly  known  as  a  Nyquist  diagram.  The  nearness  of  the 
curve  to  the  1,  -l80°  point  is  an  indication  of  the  relative  stability 
of  the  system.  The  frequency  is  also  given  for  each  point  shown  on  the 
curve.  The  predicted  curve,  shown  by  a  broken  line,  was  calculated  from 
measurements  of  the  individual  components  of  the  autopilot -aircraft  loop. 
The  solid  curve,  in  this  case,  was  calculated  from  the  closed-loop  fre¬ 
quency  response  of  the  combination  measured  in  flight. 

The  closed-loop  frequency-response  amplitude  and  phase -angle  curves 
are  shown  on  the  middle  of  each  page  plotted  against  frequency.  The 
predicted  curve  was  calculated  from  the  predicted  open-loop  response. 

The  solid  curve  was  obtained  directly  from  flight  measurements  of  the 
closed-loop  frequency  response. 

The  transient  responses  fqr  a  unit  step  input  are  shown  at  the 
bottom  of  each  page.  The  predicted  curve  was  obtained  from  the  predicted 
closed-loop  frequency  response  by  the  approximation  method  previously 
discussed.  The  solid  curve  was  obtained  from  flight  measurements  of  the 
transient  response  to  a  step  voltage  input  to  the  autopilot  equivalent 
to  l/2°  of  pitch. 

DISCUSSION 


The  three  major  aspects  of  the  analysis  to  be  discussed  are  the 
agreement  between  experimental  and  predicted  results,  the  effect  of 
displacement  and  rate  of  displacement  feedback,  and  the  effect  of 
airspeed. 

Of  primary  interest  is  the  degree  of  agreement  between  the  flight- 
measured  responses  of  the  autopilot-aircraft  combination  and  those 
predicted  from  the  component  transfer  functions.  To  serve  as  a  basis 
of  comparison,  the  important  information  which  can  be  obtained  from  the 
response  plots  can  be  summarized  as  follows:  For  the  open-loop  plot, 
the  values  of  phase  margin  and  gain  margin  specify  the  performance  to  a 
certain  degree,  the  former  usually  being  the  most  critical  with  respect 
to  stability.  On  the  closed-loop  frequency-response  plot,  the  peak- 
amplitude  ratio  and  the  frequency  at  which  it  occurs  are  significant 
points.  For  the  transient  response  there  are  several  characteristics  of 
interest,  such  as  the  response  time  (time  to  first  reach  the  desired 
value),  peak  overshoot,  cycles  to  damp  to  a  certain  fraction  of  the 
final  displacement,  period  or  frequency  of  the  oscillations,  etc. 

In  general,  it  was  found  that  agreement  between  measured  and 
predicted  responses  was  satisfactory  for  most  cases,  but  was  unsatis¬ 
factory  when  rate  signal  was  used  with  a  displacement  setting  that 
results  in  a  nearly  unstable  response  without  rate  signal. 
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The  responses  for  the  condition  of  zero  rate  and  various  sensi¬ 
tivities  and  airspeeds  are  presented  in  figures  10(a),  (b),  (c),  and  (d). 
It  can  be  seen  that  the  agreement  between  measured  and  predicted  values 
is  very  good  for  a  sensitivity  of  42  percent  at  the  three  airspeeds. 

The  comparison  at  200  knots  (fig.  10(c))  is  an  example  of  the  best 
agreement  obtained  in  the  analysis  and  is  considered  to  be  well  within 
experimental  error.  For  each  of  these  cases  at  42  percent  sensitivity, 
phase  margins  agree  within  about  5°,  and  gain  margins,  peak  amplitude 
ratios,  and  resonant  frequencies  are  within  a  few  percent.  However,  it 
may  be  noticed  in  all  these  figures  that  the  transient-response  agree¬ 
ment  is  not  as  good  as  that  for  the  frequency  responses.  The  discrep¬ 
ancy  is  actually  between  the  flight  sinusoidal  and  transient  data  since 
the  transient  peak  obtained  in  flight  is  higher  than  would  be  obtained 
from  an  analysis  of  the  corresponding  sinusoidal  response.  A  possible 
cause  of  this  difference  may  be  drift  of  the  vertical  gyro  over  the 
duration  of  the  transient  which  may  result  in  errors  when  normalizing 
the  transient-response  plot  about  the  final  steady-state  value. 

For  a  sensitivity  of  63  percent  (fig.  10(d)),  there  is  a  dis¬ 
crepancy  between  the  frequency-response  curves  which  is  greatest  at  the 
low  frequencies  on  the  open-loop  polar  diagram.  This  difference  may 
partially  be  explained  by  backlash  in  the  elevator-linkage  system. 

As  the  control  gearing  is  reduced  ( increased  servo  sensitivity) ,  the 
control-surface  deflection  becomes  smaller,  being  smallest  at  low  fre¬ 
quencies  where  the  airplane  is  closely  following  the  pitch-input  signal. 
Hence,  the  backlash  region  becomes  a  greater  percentage  of  the  total 
surface  deflection.  Thus,  the  control  gearing  would  be  effectively 
reduced  at  the  lower  frequencies  where  surface  deflection  is  smallest. 

An  increase  in  gearing  applied  to  the  flight  open-loop  response  curve 
would  bring  it  into  agreement  with  the  predicted  curve  with  the  excep¬ 
tion  that  the  frequencies  would  not  correspond. 

F°r  values  of  gearing  the  agreement  was  considered  excellent. 

Although  a  figure  is  not  presented  for  a  sensitivity  of  24  percent,  the 
predicted  ciosed-ioop  response  had  an  amplitude  ratio  peak  of  about  12 
w  ich  indicates  a  condition  dangerously  close  to  instability.  The  flight 
response  was  just  about  at  the  point  of  neutral  stability.  In  response 
to  a  transient,  the  combination  would  oscillate  sometimes  with  increas¬ 
ing  amplitude  of  oscillation  and  sometimes  with  a  very  gradual  decay  in 
amplitude  of  oscillations.  An  experimental  frequency  response  could  not 
be  obtained  because  of  the  tendency  to  break  into  instability. 

The  responses  with  rate -of -displacement  feedback  in  addition  to 
displacement  feedback  are  presented  in  figures  10(e)  through  10(i).  The 
agreement  at  a  sensitivity  of  24  percent  when  there  is  rate  feedback 
^figs.  10(e)  and  (f))  is  seen  to  be  poor;  whereas  the  agreement  at  a 
sensitivity  of  42  percent  (a  lower  gearing)  with  rate  feedback 
(.figsi  10(g),  (h),  and  (i))  is  good,  particularly  with  respect  to  phase 
margin.  As  previously  mentioned,  the  sensitivity  setting  of  24  percent 
is  one  which. results  in  a  neutrally  stable  response  without  the  addition 
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of  rate.  The  nearness  of  the  over-all  system  to  instability  at  this 

•  sensitivity  for  the  cases  with  rate  feedback  is  discussed  in  more  detail 
later  and  provides  a  clue  to  the  discrepancies.  The  system  behavior 
becomes  quite  critical  under  these  conditions  and  small  changes  in  com¬ 
ponent  values  may  result  in  large  changes  in  the  system  response.  It 
should  be  noted,  however,  that  the  curves  shown  in  figure  10(f)  repre¬ 
sent  about  the  best  transient  response  obtained  in  flight.  The  error 
voltage  was  calculated  in  each  of  the  above  cases  and  was  found  to  cause 
saturation  only  at  a  sensitivity  of  2k  percent  at  frequencies  close  to 
the  resonant  frequency. 

The  difficulty  in  determining  the  values  of  gearing  has  been 
discussed  previously  and  was  attributed  to  the  nonlinear  characteristics 
of  the  linkage  system.  The  possibilities  of  both  amplitude  and  phase  of 
elevator  deflection  being  distorted  relative  to  servo  displacement  and 
to  forcing  frequency  may  account  for  some  of  the  discrepancies  encoun¬ 
tered  in  the  predicted  results. 

A  second  major  point  of  interest  in  the  results  is  the  effect  on 
the  autopilot-aircraft  response  of  changes  in  displacement  and  rate  of 
displacement  feedback.  As  was  mentioned  previously,  a  sensitivity  of 

*  24  percent  with  zero  rate  signal  resulted  in  a  highly  oscillatory 
response  which  frequently  broke  into  instability.  Figures  10(b)  and  (d) 
show  the  effect  of  increasing  the  sensitivity  (decreasing  the  gearing) . 

As  is  to  be  expected,  the  lower  the  gearing,  the  higher  the  damping  but 
the  longer  the  response  time.  Flight  tests  bore  out  the  prediction  that 
the  response  would  be  very  slow  if  the  gearing  were  reduced  to  the  point 
of  giving  satisfactory  damping,  say  a  value  of  0.7  of  the  critical 
damping. 

The  foregoing  discussion  points  out  the  need  of  additional  stabi¬ 
lization  factors  and,  hence,  rate  of  displacement  signal  was  added  to 
the  feedback.  It  can  be  seen  from  a  consideration  of  the  transient 
responses  in  figures  10(e)  and  (f)  that  as  the  rate  signal  is  increased 
the  damping  is  improved  at  no  expense  in  response  time.  A  sensitivity 
of  2k  percent  and  a  rate  of  31  percent  (fig.  10(f))  provided  the  fastest 
response  with  the  least  amount  of  overshoot  of  any  of  the  possible  set¬ 
tings  of  the  autopilot. 

It  is  important  to  note,  however,  that  consideration  must  be  given 
other  factors  in  choosing  the  most  desirable  response  based  on  the 
criteria  given  in  the  foregoing  paragraph.  A  response  obtained  with  the 
aid  of  an  autopilot  that  is  considerably  better  than  the  response  of  the 
aircraft  alone  generally  requires  a  large  amount  of  total  control-surface 
motion.  This  may  be  objectionable  from  the  standpoint  of  servo  energy 
required  or  excessive  wear  on  the  control  system. 

An  equally  or  more  important  consideration  is  the  nearness  of  the 
system  to  instability.  When  rate  or  other  derivative  types  of  feedback 
are  used,  this  nearness  is  not  shown  by  the  closed-loop  frequency  or 
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transient  response  q/g j  since  the  output  and  input  are  not  determined 
at  the  same  point  in  the  loop.  The  relative  system  stability  is  obtained 
basically  from  an  examination  of  the  open-loop  response.  A  qualitative 
idea  can  also  be  obtained  from  the  closed-loop  response  (8+0r)/0j  (with 
reference  to  the  diagram  on  p.  }.l)  /which  may  be  rewritten  as 
0(l+PrAr/kg)/0j»  From  this  it  can  be  seen  that  the  desired  q/g j 
response  must  be  multiplied  by  the  feed-back  factor  (l+PrAr/kg)  to  check 
on  the  system  stability.  A  limit  then  exists  on  the  amount  of  rate  feed¬ 
back  that  can  be  used  to  improve  the  g/g j  response.  The  real  impor¬ 
tance  of  this  consideration  is  that  for  systems  adjusted  to  obtain  the 
optimum  G/Gi  response  small  changes  in  system  characteristics  might 
be  sufficient  in  extreme  cases  to  change  the  response  from  highly  satis¬ 
factory  to  highly  unsatisfactory. 

By  reference  to  the  experimental  results  in  figures  10(e)  and  (f) 
and  by  consideration  of  the  neutrally  stable  response  for  this  sensi¬ 
tivity  setting  (24  percent)  at  zero  rate,  it  is  seen  that  as  rate  signal 
is  progressively  increased  the  flight  transient  response  is  improved  but 
the  corresponding  open-loop  response  first  becomes  more  stable  and  then 
moves  closer  to  instability.  Further  increase  of  rate  signal  during 
flight  tests  resulted  in  actual  instability  of  the  combination  as  would 
be  expected.  This  same  trend  can  also  be  observed  in  figures  10(b) , 

(g) >  (h),  and  (i) . 

A  third  point  of  interest  in  the  results  of  the  analysis  is  the 
change  in  response  with  airspeed  for  a  given  sensitivity  setting.  This 
is  shown  in  figures  10(a),  (b),  and  (c).  Examination  of  these  responses 
indicates  that  they  are  essentially  the  same  for  all  airspeeds.  The 
phase  margins  vary  from  20°  to  30°,  gain  margins  from  2-l/2  to  5,  and 
peak  amplitude  ratios  from  2  to  2-l/2.  This  close  agreement  for  differ¬ 
ent  airspeeds  may  appear  surprising  until  it  is  recalled  that  the  gear¬ 
ing  is  unintentionally  altered  in  the  favorable  direction  by  hinge 
moment  due  to  the  elasticity  of  the  control  linkage.  It  appears  that 
the  spring  constant  is  such  that  the  change  in  gearing  compensates  for 
the  change  in  the  aircraft  response  with  airspeed. 

To  illustrate  what  would  have  happened  to  the  combination  response 
if  the  gearing  had  not  changed  with  airspeed,  consider  the  values  of  the  , 
reciprocal  gearing  l/kp  obtained  from  figure  9  for  a  sensitivity  of 
k2  percent  at  85  and  200  knots  which  are  O.92  and  1.97,  respectively. 

This  represents  a  change  of  approximately  2  to  1.  Hence,  if  the  system 
were  adjusted  to  give  the  response  at  85  knots  as  shown  in  figure  10(a), 
the  open-loop  response  for  200  knots  (fig.  10(c))  would  be  increased 
radially  by  a  factor  of  2  to  1.  The  0.8-cps  point  would  then  be  moved 
out  to  the  unity-gain  circle  resulting  in  a  phase  margin  of  only  150. 

The  gain  margin  would  be  reduced  to  about  1.6.  The  resulting  transient 
response  would  then  be  too  oscillatory  and  the  system  dangerously  close 
to  instability.  The  same  would  be  true  if  the  system  were  adjusted  at 
200  knots  and  then  flown  at  an  appreciably  higher  speed.  In  this 
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particular  installation,  the  flexibility  of  the  linkage  system  probably 
saved  the  airplane  from  destruction  in  automatically  controlled  high¬ 
speed  dives* 

CONCLUSIONS 


1-  If  the  optimum  response  (not  merely  a  conservatively  stable 
response)  for  an  autopilot-aircraft  combination  is  to  be  predicted  by 
linear  analysis  procedures,  assurance  must  be  obtained  that  the  com¬ 
ponent  performances  are  actually  linear  or  very  close  to  it*  Predicted 
responses  may  deviate  considerably  from  flight  results  due  to  a  combi¬ 
nation  of  several  small  nonlinearities  such  as  in  control  linkage  lost 
motion,  amplifier  saturation,  etc* 

2.  The  over -all -system,  open-loop  response  must  be  inspected  for 
the  nearness  of  the  system  to  instability.  Even  though  the  transient 
response  of  particular  interest  appears  well  damped,  the  system  may  be 
nearly  unstable;  thus  a  slight  change  in  component  performance  might 
result  in  instability. 

3.  The  effects  of  airspeed  on  autopilot-aircraft  stability  may  be 
compensated  for  over  the  range  in  which  the  hinge  moment  is  roughly 
proportional  to  elevator  effectiveness  by  means  of  a  simple  spring 
mechanically  linking  the  control  surface  to  the  servo  actuator.  Dynamic 
effects  of  the  spring  must,  of  course,  be  considered  in  the  design. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  Sept.  20,  1951 
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APPENDIX 


ILLUSTRATIVE  EXAMPLE  OF  CALCULATIONS 


A^/Pffd  130  knots^  sensitivity  24  percent,  and  rate  8  percent 
J.fig.  10(e))  was  selected  as  the  condition  to  illustrate  the  calcula¬ 
tions  necessary  for  a  comparison  of  predicted  and  flight  closed-loop, 
open-loop,  and  transient  responses.  This  setting  was  chosen  because, 
since  rate  feedback  is  present,  it  illustrates  in  the  most  general 
iasnion  the  calculations  required. 

The  open-loop  and  closed-loop  responses  are  calculated  in  this 
example  only  for  a  single  frequency  (f  =  0.8  cps)  since  the  procedure 
lli  be  identical  throughout  the  frequency  band.  However,  due  to  the 
nature  of  the  method,  a  complete  transient  calculation  is  carried  out. 

.  during  flight,  records  were  obtained  of  the  pitch  response  of  the 
stabilized  airplane  to  sinusoidal  and  step-inputs  of  known  magnitude. 
The  closed-loop  flight  response,  g/qt,  was  obtained  from  the  record  of 
the  attitude  of  the  airplane  G  in  response  to  the  sinusoidal  inputs: 

-  =(L  V-%  0 

VVj/  \®x J  VI 


At  f  =  0.8  cps. 


0  =  0.575  /-1970 

VJ  =  0.25  volts 

kg  =  O.51  volts/  degree 

Substituting  these  values  in  the  above  equation  gives 

—  =  1.17  /-197° 

Gj  4 - 

The  experimental  open-loop  response  AL  is  computed  from  e/er 
by  means  of  equation  (5):  0/t,I 
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where  the  feed-back  factor  ( 1+P r-^-r /kg )  can  be  evaluated  by  equation  (6) 
from  the  rate  and  no-rate  servo  responses  (fig.  7(b)),  where 


At  f  =0.8  cps. 


AL  =  Z-197°) (l.^3  / 38°) 

1-(1.17°  /-197°)(1.53 


0.66  / -1731 


_  The  Predicted  open-loop  response  was  obtained  by  multiplying  the 
dimensionless  autopilot-with-rate  response,  the  aircraft  responsl  and 
the  proper  gearing  (equation  (7)): 


AL  =  kpAp  '(.  1+ 

V  kn* 
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at  f  =  0.8  cps, 

Aq  =  O.39  /- 197°  (from  fig.  5) 

kp  =  1.52  (from  fig.  9) 

Ap  ^1+  =  1*68  /7.O0  (from  fig.  7(b)) 

Al  =  1.52  (1.68  /7.0°)(0.39  /-157°)  =  O.98  /-19QO 

The  predicted  closed- loop  response  in  terms  of  the  open-loop 
response  and  the  rate  factor',  equation  (9),  is: 


e 


al 


1+ 


PrAr 


:g 


1+Ai 


■) 


0.98  /-IgO1 


(1.53  /38°) (I+O.98  /-190o) 


=  1.24  /-114° 


The  error  voltage  of  the  autopilot-airframe  combination  was  calcu 
lated  using  the  closed-loop  no-rate  servo  response  and  the  open-loop 
flight  response  by  means  of  equation  (10). 


vec|  =  Vl  /"  !+Rg-2R  cos  €f 

V  1+1  At  l2+?l  A-r  I 


i+|al|2+2|al|  cos  eq, 

at  f  =  0.8  cps, 

R  bi  - 1,1  Z  -31°  (from  fig.  7) 

Al  =  0.67  / -167° 


vec  =  0.25  / 1+1,21  ~JLl§2  =  0. 

v  I+O.45  -  I.31 
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The  predicted  transient  response  was  obtained  by  FloydTs  Method 
which  is  outlined  on  page  14.  For  airspeed  130  toots ,  rate  8  percent, 
sensitivity  24  percent,  Re  H(  jw),  which  in  this  case  is  Re  d/0 i(jco), 
was  plotted  as  a  function  of  frequency  f  and  the  resultant  curve  was 
approximated  by  straight-line  segments  as  shown  in  figure  11.  The 
values  of  r,  f,  and  the  corresponding  00  are  tabulated  as  follows: 


ro  = 

1.00 

^a 

rl  = 

.75 

*b 

r2  = 

.35 

fC 

r3  = 

.60 

fd 

r4  “ 

.95 

fe 

ff 

0.34 

a;  a  =  2.14 

•  53 

coij  =  3*34 

.64 

coc  =  4.03 

.83 

ud  =  5.22 

1.004 

(0  e  =  6.30 

130 

Mf  =  8.17 

1.75 

co  g  =  10  •  10 

Since  all  of  the  component  trapezoids  must  be  measured  from  the 
origin.  Re  H(jco)  was  approximated  by  the  sum  of  the  trapezoids  (with 
proper  sign  affixed)  as  shown  in  figure  12.  Upon  evaluation 


7 


n=  1 


yielded  the  impulse  response  shown 
yielded  the  step  response  shown  in 


in  figure  13,  which,  upon  integration, 
figure  10(e). 
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Figure  1.—  Photograph  of  the  airplane. 
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Aircraft  attitude, 


fograph  records  obtained  during  longitudinal  dynamic 
tests  of  the  stabilized  test  airplane . 


Amplitude  ratio,  /Aff/ 
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Frequency,  cps 
(a)  Amplitude  ratio . 


Figure  5-  Airplane  frequency  response  in  pitch . 


Phase  angle, 
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(b)  Phase  angle. 


Figure  5-  Concluded . 
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(a)  Rate  O,  various  sensitivities,  ±0.U5  volt  (&  0.25°)  input . 
Figure  7.-  Autopilot  response. 
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(b)  Sensitivity  24  percent,  various  rates,  ±  0.25 0  input. 


Figure  7.-  Concluded . 
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Figure  9- Reciprocal  gearing,  l/kp,  versus  hinge  moment  for  various 

sensitivities . 
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Figure  tO  -  Continued . 
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Figure  //.-  Approximation  of  real  part  of  dosed -loop  response,  Re  -f^  (ju),  by  several  trapezoids 

for  130  knots,  sensitivity  24  percent,  and  rate  8  percent. 


